Reentrant spin-glass transition in a dilute magnet 



S. Abiko, S. Niidera, and F. Matsubara 

Department of Applied Physics, 
Tohoku University, Sendai 980-8579, Japan 

(Dated: February 2, 2008) 

We have performed a large scale Monte Carlo simulation of a dilute classical Heisenberg model 
with ferromagnetic nearest neighbor and antiferromagnetic next-nearest neighbor interactions. We 
found that the model reproduces a reentrant spin-glass transition. That is, as the temperature is 
decreased, the magnetization increases rapidly below a certain temperature, reaches a maximum 
value, then ceases at some lower temperature. The low temperature phase was suggested to be a 
spin-glass phase that is characterized by ferromagnetic clusters. 
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Reentrant spin-glass (RSG) transition is a well-known 
phenomenon of spin glasses (SGs). The RSG transition 
is found near the phase boundary between the SG phase 
and the ferromagnetic phaseO, @|- As the temperature 
decreases from a higher temperature, magnetization in- 
creases, then ceases at a lower temperature. Finally, the 
SG phase is realized. The phenomenon was first con- 
sidered as a phase transition between a 'ferromagnetic 
phase' and a 'SG phase'|3|. However, neutron diffraction 
studies have revealed that the 'SG phase' is characterized 
by ferromagnetic clusters 0, 0, 0. Now the RSG transi- 
tion is believed to be a re-entry from a ferromagnetic 
phase to a frozen state with ferromagnetic clusters. 

The mechanism responsible for this reentrant transi- 
tion has not yet been resolved. Two ideas have been 
proposed for describing the RSG transition: (i) an 
infinite-range Ising bond model, and (ii) a phcnomeno- 
logical random field concept. Sherrington and Kirk- 
patrick solved the infinite-range Ising bond model using 
a replica technique 0. They predicted the occurrence of 
the RSG phase transition before the RSG transition was 
observed experimentally 0. However, successive stud- 
ies of the model revealed that the RSG phase transition 
does not occur Q, even when the vector spins substitute 
for the Ising spins|||. On the other hand, the random 
field idea was proposed to explain experimental obser- 
vations of neutron scattering functions |J, !! G3. The 
essential point of that conception is that the system is 
decomposed into a ferromagnetic part (FM part) and a 
part with frustrated spins (SG part). At low tempera- 
tures, the spins of the SG part will yield random effec- 
tive fields to the spins of the FM part. The FM order, 
which grows at higher temperatures, vanishes because 
of a random field effect [Tl| . Later, it was pointed out 
that not random field but random anisotropy brings a 
RSG-like phenomenon [H El • 

In a FM model with ran- 
dom anisotropy, as the temperature is decreased from 
a high temperature, there first appears a correlated SG 
(CSG) phase that is characterized by a well developed 
FM spin correlation. Then the CSG phase will become 
transformed into a speromagnet, i.e., a random SG-like 
phase, just like the RSG transitionjlj- Nevertheless, no 



theoretical evidence has yet been presented on this idea 
in a microscopic point of view. In the last two decades, 
computer simulations have been performed extensively 
to solve the RSG transition in various models such as 
short-range bond models Hi EU, short-range site 
model s Il7l Ha and a Ruderman-Kittel-Kasuya-Yoshida 
modelpjj, |2(j. However, no realistic model has been 
shown to reproduce vanishing magnetization that grows 
at higher temperatures. 

This Letter reports a dilute classical Heisenberg model 
that reproduces a RSG transition. That is, the mag- 
netization increases rapidly below a certain tempera- 
ture, then disappears at some lower temperature. The 
low temperature phase was suggested to be a SG 
phase, which is characterized by ferromagnetic clus- 
ters. It should be emphasized that no novel mecha- 
nism is necessary to reproduce the RSG transition. The 
model studied here is a natural one that was proposed 
experiment allyj2l| and investigated using a Monte Carlo 
(MC) method [l7| . A large-scale MC simulation revealed 
the nature of the model. We believe that the system size 
is crucial for vanishing the magnetization, because the 
SG phase at low temperatures is composed of ferromag- 
netic clusters. 

We start with a dilute Heisenberg model with compet- 
ing nearest and next-nearest neighbor exchange interac- 
tions described by the Hamiltonian: 

nn nnn 

H= - }^ JiXjXjSj ■ Sj + y] J 2 XfcX;gfc • Si, (1) 

where Si is the classical Heisenberg spin of \Si\ = 1; 
Ji(> 0) and J 2 (> 0) respectively represent the nearest 
neighbor and the next-nearest neighbor exchange inter- 
actions; Xi — 1 or when the lattice site i is occupied 
respectively by a magnetic or non-magnetic atom. The 
average number of x(= (xi)) is the concentration of a 
magnetic atom. Note that an experimental realization 
of this model is Eu^Sri_^S|2l||. in which magnetic atoms 
(Eu) are located on the fee lattice sites 22). Here, for sim- 
plicity, we consider the model on a simple cubic lattice 
with J 2 = 0.2 Jn|2^. 



2 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

T/Jj 

FIG. 1: Magnetizations M in the 32 x 32 x 32 lattice for 
various spin concentrations x. 



A computer simulation was performed using a conven- 
tional heat bath MC method. The system was cooled 
gradually from a high temperature (cooling simulation). 
We calculated the magnetization M defined as M = 
[(M(s))], where M(s)(= | XiSi\) is the magnetization 
at the sth MC step, and (■ • ■ ) represents a MC aver- 
age and [• • • ] a sample average. Here, for larger lattices, 
200000 MC steps (MCS) were allowed for relaxation; data 
of successive 200000 MCS were used to calculate average 
values. We will show later that these MCS are suffi- 
cient for studying equilibrium properties of the model at 
a temperature range within which the RSG behavior is 
found. We made a simulation for x > 0.70. We treated 
lattices of Lx Lx L (L = 8 — 48) with periodic boundary 
conditions. Numbers N s of samples with different spin 
distributions arc N s = 1000 for L < 12, N s — 600 for 
L = 16, N s = 200 for L = 24 and 32, and yV s = 80 for 
L — 48. We measured the temperature in the unit of J\ 

(*B = I)- 

Figure 1 shows temperature dependencies of the mag- 
netization per spin M/xN for various x, where N(= L 3 ) 
is the number of lattice sites. For x = 1, as the tempera- 
ture decreases, M increases rapidly below a temperature 
revealing the occurrence of a ferromagnetic phase. As 
x decreases, M exhibits an interesting behavior. In the 
range of 0.78 < x < 0.85, M once increases, reaches a 
maximum value, then decreases. Such behavior of M is 
reminiscent of the occurrence of the RSG transition. 

To examine this phenomenon, we made detailed stud- 
ies in the case of x = 0.80. First we note that we 
performed a complementary simulation. That is, start- 
ing with a random spin configuration at a low tempera- 
ture, the system is heated gradually (heating simulation) . 
In addition, we investigated M in the ground state for 
smaller lattices (L < 16) having used a hybrid genetic 
algorithm|24j. Figure 2 shows temperature dependencies 
of M in both cooling and heating simulations for vari- 
ous L together with the value of M in the ground state. 
For T > 0.1 Ji, data of the two simulations almost coin- 
cide mutually, even for large L. We thereby infer that M 
for T > O.Ui are of thermal equilibrium. On the other 
hand, for T < 0.1 Ji, a great difference in M is seen be- 
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FIG. 2: Magnetizations M for x = 0.80 in the L x L x L 
lattice. Open symbols indicate M in the cooling simulation 
and filled symbols that in the heating simulation. Data at 
T — indicate those in the ground state. 
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FIG. 3: Binder parameters §l for x = 0.80. 



tween the two simulations; estimation of the equilibrium 
value is difficult. We speculate that the heating simu- 
lation gives a value of M that is similar to that in the 
equilibrium state because the data in the heating simu- 
lation seem to connect to those in the ground state. 

Figure 2 shows that the lattice size dependence of M is 
remarkable. For smaller L, as the temperature decreases, 
M decreases slightly at very low temperatures. The de- 
crease is enhanced as L increases. Consequently, a strong 
size-dependence of M is indicated for T < O.Ui. Particu- 
larly in the ground state, M apparently decreases rapidly 
as L increases. These facts imply that M disappears at 
low temperatures as well as at high temperatures. 

We investigated the Binder parameter g^, 



9l = (5 



JW0 4 )] wo 
[<M( S ) 2 >] 2i/ ' 



(2) 



to examine whether or not the ferromagnetic phase ac- 
tually occurs within an intermediate temperature range. 
Figure 3 shows g^'s for various L|25j. We can see that 
<7l's for different L cross at two temperatures Tc and Tr 
(< T c ). The cross at T c /Ji = 0.265±0.010 is a usual one 
that is found in the ferromagnetic phase transition|26j. 
That is, for T > Tc, gh for a larger size is smaller than 
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FIG. 4: The square of the magnetization [(M(s) 2 )]/xN vs. L. 
The straight line at T = O.Ui is obtained in the least square 
fitting for smaller lattices of L < 16. For T = 0.07Ji, open 
symbols indicate those in the cooling simulation and filled 
symbols represent those in the heating simulation. Data in 
both simulations deviate downward from the straight line. 



that for a smaller size; for T < Tc, this size dependence 
in gL is reversed. On the other hand, the cross at Tr is 
strange. For T < Tr, for a larger size again becomes 
smaller than that for a smaller size. Interestingly, the 
cross for different gL occurs at almost the same temper- 
ature of T R /Ji = 0.125 ± 0.005(2^. These facts suggest 
that, as the temperature is decreased beyond Tr,, the 
ferromagnetic phase, which occurs below Tc, disappears. 

Does the spin correlation function below Tr decay 
rapidly as the distance between two spins increases? If it 
decays according to the power law, then [(M(s) 2 )]/xN oc 
i 2 ~ r ' with rj being the decay exponent of the spin corre- 
lation function. We plot, in Fig. 4, {(M{s) 2 )]/xN as a 
function of L in a log- log form at various temperatures. 
In fact, at a temperature near below Tr,, data seem to 
lie on a straight line; the spin correlation function will 
decay algebraically. Different L dependences are found 
above and below this temperature. At higher tempera- 
tures (T > Tr), as L is increased, data deviate upward 
from the straight line being compatible with the fact that 
the ferromagnetic long range order occurs at these tem- 
peratures. On the other hand, data deviate downward at 
lower temperatures. That is, the spin correlation func- 
tion will decay exponentially, similarly to that in the SG 
phase. 

Is the SG phase realized at low temperatures? A con- 
vincing way of examining the SG phase transition is a 
finite size scaling analysis of the correlation length, 
in samples of different sizes L [2]], [2^] • Data for the dimen- 
sionless ratio £,l/L are expected to intersect at T = Tsg- 
Here we consider the correlation length of the SG com- 
ponent of the spin, i.e., Si(= Si — m) with m being the 
ferromagnetic component of m = '^2, i XiSi/{xN). We 
performed a cooling simulation of a two-replica system 
with {Si} and {Ti} 29] . The SG order parameter, gen- 
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FIG. 5: The SG correlation length £l divided by L. The data 
intersect at T ~ O.lOJi, implying that there is a SG transition 
at this temperature. The inset shows a typical example of the 
scaling plot for L > 12. 

eralized to wave vector fc, q >J ' v (k), is defined as 

^(fcH^E^Te^, (3) 

i 

where fj,,u = x,y,z. From this, we determine the wave 
vector dependent SG susceptibility xsg(&) by 



XSG 



(k) = xNj2i(W(k)\ 2 )} 



ft.iy 



The SG correlation length is then determined from 



i 



Xsg(0) 



2 sin(fc min /2) 



1) 1/2 , 



(4) 



(5) 



where fe min = (2tt/£,0, 0). 

Figure 5 shows the temperature dependence of ^l/T 
for various L. In fact, £l/T for different L intersect at 
T ~ O.lOJi. In particular, data for L > 12 are scalable 
on the assumption that Tsg/<T = 0.105 ± 0.003. On 
the basis of this fact together with the rapid decay of 
the spin correlation function, we inferred that the SG 
phase is realized at low temperatures. The SG transition 
temperature Tgc estimated here is slightly lower than 
Tr 30]. However, the possibility of Tr = Tsg cannot 
be ruled out, because the treated lattices of L < 20 for 
estimating Tsg are not large enough. 

We considered the spin structure. Figures 6(a) and 
6(b) show typical results for it. In the ferromagnetic 
phase (Fig. 6(a)), although the spin arrangement is con- 
siderably modulated, a ferromagnetic spin correlation ex- 
tends over the lattice. On the other hand, in the SG 
phase (Fig. 6(b)), we can see that the system breaks 
up to yield ferromagnetic clusters with a linear size of 
l c ~ 7 31]. This result is compatible with the size de- 
pendence of M shown in Fig. 2. In the small lattice 
with L = 8, M does not exhibit a marked decrease at 
T < Tsg- The decrease becomes drastic when L is in- 
creased. This fact is further evidence that the system is 
divided into ferromagnetic clusters. We suggest that the 
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SG phase of the model is characterized by ferromagnetic 
clusters. This concept is compatible with experimental 
observations 0, . 

In summary, we found a model that settles the most 
important issue of the RSG transition. Other important 
issues remain unresolved. Why does the ferromagnetic 
phase disappear at low temperatures? Does this dilute 
model exhibit the same behavior as that found in the 
bond SG model[2^, |32? Does the chiral glass phase 
transition simultaneously occur at Tsg|11|? We intend 
the present model as one means to solve those and other 
remaining problems. 
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(a) T = 0.19Ji 



(b) T = 0.04Ji 



xxx* xx\-,\xxxx 

* x. XXX XXXXXX 
X*. XXXXX* XXXXX* 

XXX xxxxxxxxxx ^ 

*; \ x x x xxx xxx t 

* t X-.XXX x * 
1 * t t ■ XXXXXXX * 

* M X ■ I 

< 1 XXXXK-.X \t t 

WW xxxx\\xxx t 
t t WWWW WW 
w w.\w\w\w.t 1 1 

W t »; t H X X << s <t »< •> 

^ . XXXX a .H- XXXX 

»x*,xxxx im^Ml \ ^ 

^ », x ^^M*-. xx 
x \xxx * ■ I 

* Hvxxx 



1 1 t t k\xxxxx > t t 

WW Wxxxxx t f t 

•••..-». » » • »a 

1 M * H xxxn < * I 

hhi i », t 

\ t * xxx *^>i^ 

1 ft xxxx 1 », n \ 
». i t ^\\w n t 
ttttt xxxxx>s *, \ 
f 1 1 1 1 \ xxxxxx\x 1 

tttt \ XXXXX-- 
t t x XXX"*- 

HI' - ----- » x x 1 

t \ t - 

i, ^ - - » 

- - - - - - - \ ^ * 

XXXXxxx-ox---. 

X XXXXX-. -.XX X *\ 




- • «--x--XX U\ 

""Wi ' ' «-x\\M 
m H ■«•• • • » >'-<>.v^HM 

• tttttt -.-■.MttHU'-M 
«- tttttt -,xxx uttn It 

• _»........ ■ 

• • -• »<>* ft • .« • 

»■■_ ft'* 5 tttt 
•H ■ ' '- «-'-( ft t 
*^Hx f f t x - - - - - 

t * / MltHH* •»■*-'■ 
1 (HUHI *< W 

• $ * f t ^ f 






FIG. 6: Spin structures of the model for x = 0.80 on a plane of the 32 x 32 x 32 lattice in (a) the ferromagnetic phase and (b) 
the spin glass phase. Spins represented here are those averaged over 10000 MCS. The positions of the non-magnetic atoms are 
represented in the white. 



